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This paper proposes a multi-choice goal programming model for dealing with the capacity expansion 
planning problem of the renewable energy industry. This model involves decisions regarding the optimal 
mix of different plant types, location selection and other criteria. Different types of plants should be 
located in appropriate places so as to minimize the total deviations from predefined goals concerning 
power generated, investment cost, emission avoided, jobs created, operation and maintenance costs, 
distance security, and social acceptance. The proposed method is superior to the goal programming 
model proposed by Ramón and Cristóbal, in that it can avoid underestimation of aspiration level, expand 
the potential feasible region, and achieve findings more closely approach actual conditions. In addition, 
the social acceptance of the renewable energy planning problem in Taiwan is modeled by the MCGP to 
demonstrate its usefulness. 
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1. Introduction 


Because of continuous economic growth and industrialization 
worldwide and the soaring energy demands in the five major 
developing countries of Asia (China, India, Indonesia, Korea and 
Thailand), the demands on fossil energy for industry, transporta- 
tion and daily life are becoming increasingly large. Table 1 shows 
the global reserves and the projected number of years of supply 
remaining for major fossil fuels [1,2]. Even with recent discoveries 
of new reserves, such as shale gas, tight sands and shale oil in 
North America, as well as shale oil in many other areas, fossil fuels 
are gradually being depleted. International oil and coal prices 
exceeded US $102 per barrel and US $75 per ton in April 2010 [3,4]. 
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Since Taiwan has almost no energy resources, 97.92% of the total 
energy requirement depended on imports in 2013 and almost all of 
the fossil fuel is imported from politically unstable areas, such as 
the Middle East, so Taiwan's energy security has become increas- 
ingly vulnerable. Energy imports also increased from 3.88% of 
Taiwan's GDP in 2002 to 14.55% in 2012, which creates a large 
drain on Taiwan's economy. In addition, traditional uses of fossil 
fuels aggravate environmental pollution and the greenhouse effect. 
Thus, the development of renewable energy is an increasingly 
important issue for both developing and developed countries. 
Taiwan's energy supply and consumption in 2012 are listed in 
Fig. 1. Compared with 2011, petroleum is 3.4% higher, renewables 
are 0.1% higher, nuclear is 0.2% higher, natural gas is 0.7% higher, 
and coal is 4% lower. This means that Taiwan is gradually improving 
with the use of cleaner fuels to reduce carbon emission from power 
generation. Since the 2011 nuclear disaster at Fukushima, the 
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Table 1 
Global reserves and availability of major energy resources. 


Category item Oil (+ oil sanss) 
Total reserves (end of 2011) 
Yield (2011) 

Available years 


1481 billion barrels 
30 billion barrels 
49 years 


2012 Energy Supply and Consumption 


Petroleum 48.0% gaze 


Renewables 1.996 


Natural Gas 12.1% 
Nuclear 8.335 — ————————— 


Transportation 11.996 
Service 11.0% 


Industrial 38.2% 
Residential 10.9% 


Energy 7.1% 


Non-energy Use" 20.0% Agricultural 0.9% 


Fig. 1. 2012 Energy supply and consumption in Taiwan. 


government of Taiwan has considered a nuclear energy reduction 
policy. Overall, it is a top priority of Taiwan to develop clean, 
sustainable, and independent energy and achieve a balance between 
energy security, environmental protection, industrial competitive- 
ness, and economic development. This is particularly important for 
Taiwan's environment, with the sunny weather in southern Taiwan 
and abundant wind resources on the west coast of Taiwan. These 
issues lead to both crisis and opportunity for Taiwan in becoming a 
sustainable community. In order to encourage the development of 
renewable energy, the National Science Committee (NSC) initiated 
the flagship project “New Projects of Energy Technology” in 2008. 
Subsequently, the Ministry of Economic Affairs (MOEA) promulgated 
the “Renewable Energy Development Act (REDA)” in 2009, the 
“Incentive Program of Offshore Wind Power Demonstration System” 
in 2012, the “Million Solar Rooftop Photovoltaic Program” in 2013, 
and the “Small Wind Turbine Generation System Demonstration 
Incentives” in 2014. Many kinds of financial compensation are offered 
by the central and local governments. Further programs are also 
expected under the Bureau of Energy (BOE) for promoting renewable 
energy usage: (i) Total renewable energy capacity of 9952 MW in 
2025 (i.e., 15.5% of the total power generation installation capacity in 
Taiwan); (ii) 1000 on- and off-shore wind turbines, with the wind 
power capacity of 4200 MW in 2030 and progress on the construc- 
tion of the wind power system as shown in Table 2; and (iii) the 
"Million Solar Rooftop Photovoltaic Program" project, with photo- 
voltaic power capacity of 3100 MW in 2030. In particular, the 
Industrial Technology Research Institute (ITRI), Taiwan Power Com- 
pany, Chinese Petroleum Corporation, and others have begun to carry 
out research on renewable energy technology. From the promulga- 
tion of the REDA in 2009 to the end of 2012, the capacity of 
photovoltaic power increased from 11 MW to 222.4 MW, the capa- 
city of wind power increased from 436 MW to 621 MW. 


Natural gas Coal 


860 billion tons 
6.7 billion tons 
112 years 


196 trillion cubic meters 
3.2 trillion cubic meters 
61 years 


However, these infrastructure projects are often contested and 
face public resistance, in particular against many of the on-shore 
wind power projects in Taiwan. Therefore, public resistance should 
be considered as an important factor in a model for the optimal 
location of renewable energy facilities. The capacity expansion 
planning problem of the renewable energy industry should not 
only to determine the best renewable energy project among 
different alternatives, but also develop expansion plans for energy 
generation, transmission and distribution systems with considera- 
tion of the potential public resistance to these projects. This 
problem involves important decisions regarding the optimal solu- 
tions of different plant types, locations, and capacities. A variety of 
approaches have been employed to model the generation and 
distribution of electricity, such as of multi-criteria decision-mak- 
ing models [5-7], optimization models with geographic informa- 
tion systems for regional renewable energy development [8], and 
goal programming (GP) model for the optimal mix and location of 
renewable energy plants [9]. As a result, multi-criteria decision- 
making methods such as the Analytical Hierarchy Process (AHP), 
Elimination and Choice Expressing Reality (ELECTRE), Preference 
Ranking Organization Method for Enrichment of Evaluations 
(PROMETHEE), Technique for Order Preference by Similarity to 
an Ideal Solution (TOPSIS), and VlISekriterijumsko KOmpromisno 
Rangiranje (VIKOR) have been extensively adopted for the selec- 
tion of renewable energy problems such as wind-farm projects, 
geothermal projects and hydro-site selection | 10-16]. 

Though these methods have been widely applied to the renew- 
able energy industry, the multi-choice goal programming (MCGP) 
method [17-19] has rarely been used for the location selection 
problem of renewable energy plants and capacity expression 
decisions. MCGP is a multi-choice aspiration-level approach which 
has been applied to many real-world problems, such as housing 
selection [20], low-cost carriers’ networks [21], supplier selection 
[22,23], forestry plans [24], forestry management [25] and a 
piloted quality-management system [26]. Moreover, using MCGP 
to deal with the optimal location of renewable energy plants, can 
provide the following advantages: 


* MCGP expands the original feasible region (it may be obtained 
using Weight GP, MINMAX GP, or Lexicographic GP) to the 
potential feasible region for obtaining the appropriate solution 
for decision-makers (DMs) responsible for the selection of the 
renewable energy plant location [17]. 

* It not only can deal with quantitative issues (e.g., minimization 
of operation and maintenance costs) but also qualitative 
problems (e.g., maximization of user satisfaction and social 
acceptance) [17]. 

* [t allows DMs to set multiple aspiration levels to avoid under- 
estimating settings determined by conservative chief-executive 
officers. For example, the lower power capacity of 
113,609,530 W is obtained in the model of Ramón and Cristóbal 
(RC) [9] because the power generation goal (110 x 109 W) is 
underestimated in their model. In contrast, a higher power 
capacity of 116,587,600 W is obtained by the MCGP model 
because MCGP can reflect higher aspirations. That is, by 
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avoiding underestimation of the aspiration level, it is possible 
to have better solution in MCGP [19]. 


This study uses MCGP to deal with the selection of renewable 
energy plants, locations and the capacity expression problem by 
taking into account the above advantages and striving to achieve 
more reality-based findings. The superiority of the proposed 
method can be seen in the comparison of the solutions between 
Ramón and Cristóbal [9] and the proposed MCGP method. In 
addition, the MCGP is used to solve a wind power plant location 
selection problem in Taiwan taking into account public resistance. 


2. Ramón and Cristóbal model (RC model) 
In order to locate five renewable energy plants for power 
generation in five places located in the region of Cantabria in the 


north of Spain, Ramón and Cristóbal [9] proposed a novel GP 


Table 2 
1000 wind turbines project in Taiwan. 


Years 2012 2015 
On-shore wind 621 866 

(MW) (314 turbines) (350 turbines) 
Off-shore wind 0 15 

(MW) (4 turbines) 


Total 621 881 
(314 turbines) (354 turbines) 


Table 3 
Alternatives for electric generation. 


model to address the problem of how to locate seven candidate 
plants (alternatives) in seven locations: La Braguía, La Vega, 
Estacas, Lunada, Potes, Reinosa and Santander. The attributes 
considered when evaluating these renewable energy systems in 
this model are: power generated (W), investment cost (I), tons of 
emissions of CO; avoided per year (TCOs5/y), jobs created (J), 
operation and maintenance costs (OM), distance between plants 
(D) and social acceptance (S). Table 3 shows the data for the first 
five attributes. Table 4 shows the distance between the locations. 
Table 5 shows the social acceptance by the local population for the 
hypothesized realization of the renewable energy projects. Social 
acceptability is expressed using a scale of 1 (low acceptance) to 10 
(high acceptance). As different types of plants can be placed in 
each location, the goal is to ensure the minimum total deviation 
from the goals. These goals are given as 


- (G1) The power generated must be higher than 110 x 109 W. 
- (G2) The investment cost must be limited to 100 x 10°€/year. 


2020 2025 2030 
1200 1200 1200 
(450 turbines) (450 turbines) (450 turbines) 
600 1800 3000 
(120 turbines) (360 turbines) (600 turbines) 
1800 3000 4200 


(570 turbines) (810 turbines) (1025 turbines) 


Alternative Ww I TCO;/y ] OM 
Wind power (10 x P x 50 MW) 60,959,000 23,425,000 9,649,680 15 37,750 
Hydroelectric (P < 10 MW) 3,940,100 7,500,000 472,812 8 7250 
Hydroelectric (10 £ P £ 25 MW) 962,000 14,000,000 255,490 8 14,000 
Hydroelectric (25 < P < 50 MW) 412,000 21,035,000 255,490 12 21,000 
Solar (P > 10 MW) 306,031 5,320,000 482,856 10 1792 
Biomass (P < 5 MW) 13,612,500 9,015,000 2,524,643 15 27,100 
Biomass (P > 50 MW) 37,770,000 47,958,400 4,839,548 20 75,000 
Table 4 
Distance between locations. 
Alternative La Braguia La Vega Estacas Lunada Potes Reinosa Santander 
La Braguia - 6.5 15 39 126 61 53 
La Vega - 32 121 54 65 
Estacas 24 128 62 72 
Lunada - 143 78 56 
Potes - 130 107 
Reinosa - 75 
Santander - 
Table 5 
Social acceptance (1-10) of each plant in each location. 
Alternative La Braguía (8) La Vega (9) Estacas (10) Lunada (11) Potes (12) Reinosa (13) Santander (14) 
(1) Wind power (10 < P < 50 MW) 7 9 7 9 
(2) Hydroelectric (P £ 10 MW) 5 6 5 7 
(3) Hydroelectric (10 £ P £ 25 MW) 6 7 6 7 
(4) Hydroelectric (25 « P £ 50 MW) 4 5 6 
(5) Solar (P > 10 MW) 8 7 9 
(6) Biomass (P < 5 MW) 9 8 8 9 
(7) Biomass (P 2 50 MW) 9 8 9 
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- (G3) The emissions avoided must be higher than 18 x 10° t. 

- (G4) The jobs created must be higher than 70. 

- (G5) The operation and maintenance costs must be limited to 
350,000€J year. 

- (G6) The distance between plants must be maximized to 
2910 Km. 

- (G7) The social acceptance must be as close as possible to the 
highest level of 50. 


Fig. 2 shows the multi-source multi-sink network involving 
seven alternatives and seven locations, a flow of k units from the 
source node through alternatives and place nodes to the sink node 
represents an assignment of k plants to places in which k 
compatible couples are created. This behavior can ensure that 
alternatives and place nodes are represented by binary variables. 
Then Ramón and Cristóbal [9] proposed a GP model (see [9] for 
details) for the optimal mix and location of renewable energy 
plants. In their model, the optimal mix of different plant types, 
locations, and capacity expansion can be easily determined using 
the multi-source multi-link network analysis mechanism. The 
compromise solutions of RC model thus obtained are shown in 
Table 6. 


3. Improvement of application 


The GP approach is a well-known and popular multi-objective 
mathematical programming model. It can be accomplished with 
various types of methods such as Lexicographic GP, Weight GP, 
MINMAX (Chebyshev) GP, and Extended GP [27]. Its inherent 
flexibility allows users to take into account several conflicting 
objectives, multiple criteria, and incomplete information at same 


Places 


Alternatives 


m ed 


o 
OS" 


o 
Sr 


"A 


SY 


Fig. 2. Multi-source multi-sink network. 


Table 6 
Solutions comparison of various models. 


time for choosing the most satisfactory solution within a feasible 
region. GP was first introduced by Charnes et al. (1955) and further 
developed by Charnes and Cooper (1961), Lee (1972), Ignizio 
(1985), Tamiz et al. (1998) and Romero (2001), among others 
[17]. In practice, the GP model minimizes Y;(d; Hd; ), subject to 
constraints f;x)—d; +d; =g; where xeX CR", d; and d; are 
positive and negative deviations of the achievement function f;(x) 
from the scalar aspiration level g;. In other words, the purpose of 
GP is to minimize the deviations between the achievement of 
goals and the scalar aspiration levels (it may be predefined by DMs 
for their problems). 

In contrast, the role of MCGP is to aid a DM in expanding the 
original feasible region, which may be obtained by GP, to the 
potential feasible region to find more and better solutions. 
The MCGP model minimizes Yd; +d; te; --ej ), subject to 
constraints fix% -df +d; =y; Yi-ej ej =8i max OT Zi, mim 
where xeX CR", d and d; are positive and negative deviations 
of the achievement function f;(x) from the vector aspiration level 
Yi, ej and e7 are positive and negative deviations of the vector 
achievement level y; from the scalar aspiration level g; max or 
Zi, min- Here, gi max (Zi, min) is the sufficiently large value for the case 
of the more (less) the better in real problems. In order to more 
closely approach reality with a potential solution than the com- 
promise solution obtained by Ramón and Cristóbal [9], the 
aspiration for each goal in the model of Ramón and Cristóbal 
should be changed as follows: 


- (G1) Power generated must be higher than 110 x 109 W and 
the higher the better. 

- (G2) Investment cost must be limited to 100 x 10°€/year and 
the lower the better. 

- (G3) Emissions avoided must be higher than 18 x 106 t and the 
higher the better. 

- (G4) Jobs created must be higher than 70 and the higher the 
better. 

- (G5) Operation and maintenance costs must be limited to 
350,000€/year and the lower the better. 

- (G6) Distance between plants must be maximized to 2910 Km 
and the greater the better. 

- (G7) Social acceptance must be as close as possible to the 
highest level of 50 and the higher the better. 


Table 7 

Places to locate the plants of various models. 
Alternative MCGP-MSMS model 

1) Wind power (10 < P < 50 MW) 


2) Hydroelectric (P £ 10 MW) 


( Estacas (10) 
( 

(5) Solar (P 2 10 MW) 

( 

( 


Reinosa (13) 
Potes (12) 
Lunada (11) 
Santander (14) 


6) Biomass (P < 5 MW) 
7) Biomass (P 2 50 MW) 


W (G1) 1(G2) TCOsJy (G3) 
Goals 110 x 10° 100 x 10° 18 x 10° 
RC model 113,609,530 99,718,400 17,752,217 
MCGP-MSMS mdel 116,587,600 93,218,400 17,969,540 
MCGP-MSMS-normalized 116,587,600 93,218,400 17,969,540 
MCGP-MSMS-normalized? 113,059,500 106,753,340 17,752,220 
MCGP-MSMS-normalized^ 116,587,600 93,218,400 17,969,540 


* 10 times as a weight add to job created goal. 
P 10 times as a weight add to social acceptance goal. 


J (G4) OM (G5) Distance (G6) Social acceptance (G7) 
70 350,000 2910 50 
68 155,642 2323 43 
68 148,892 1043 42 
68 148,892 1043 42 
72 162,642 1043 41 
68 148,892 948 43 
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For treating these revised goals with new aspiration levels, this 
study proposes a novel method, referred to as the MCGP with a 
multi-source and multi-sink network (MCGP-MSMS) model for 
solving the capacity expansion planning problem of the renewable 
energy industry (see the Appendix for details). This MCGP-MSMS 
model is then solved using LINGO [28] to obtain the solution, as 
shown in Tables 6 and 7. Table 6 shows the deviations from the 
goals, while Table 7 and Fig. 3 show the selected locations for each 
plant. As seen in Table 6, goals G1, G2 and G5 are achieved, while 
goals G3, G4, G6 and G7 are not in both RC and MCGP-MSMS 
models. 

It is noteworthy that the power generated (G1), investment 
cost (G2), emissions avoided (G3), and operation and maintenance 
cost (G5) achieved in MCGP-MSMS model are better than the 
results obtained by the model of RC [9]. Compared with the RC 
model, the power generated is 2.62% higher, investment cost is 
6.97% lower, emissions avoided are 1.22% lower, and operation and 
maintenance costs are 4.53% lower in the MCGP-MSMS model. 
Obviously, the underestimation problem can be overcome by the 
proposed MCGP-MSMS model and the original feasible region in 
the GP has been expended to the larger potential feasible region in 
MCGP-MSMS [17]. The features of the MCGP-MSMS and GP are 
shown in Fig. 4. 


Alternatives 


Fig. 3. The solution of MCGP-MSMS model. 


— (C) Ideal solution 


; MCGP-MSMS 
< solution 


Lexicographic 
GP solution 


Potential 


feasible 


Weight GP 
solution 


Original 
feasible 


region 


Fig. 4. Possible feasible regions in various GP methods. 


Deviational variables d;", dj , e;* and e; in the model of RC [9] 
or MCGP-MSMS can vary in different ranges, resulting in an 
unintentional bias towards objectives with a large magnitude. 
For example, the goals of power generated and investment cost in 
the model of RC [9] require a relatively large value to obtain a large 
value in the objective function. This forces these two goals to be 
given higher priority than other goals to reach the target in the 
optimization model. Therefore, the use of a normalization techni- 
que to overcome this problem is an important issue in multiple- 
criteria decision-making. In order to solve this problem, the utility 
normalization technique (Chang [17]) is slightly modified as 
follows: 

Let n, =[0,U,] and p, 2[0, Vj] Yk where U; and V, are the 
upper bounds of n; and pẹ, respectively. D, = Max (U,, Vk} Y k and 
Lj - Min (U,, Vk} Y k. The normalized weights are easily calcu- 
lated as w,— (L;/Lj 4- Dy), Wj = (Dy /L; +D,), and w;= (wjU,/Ui) 
(ie {1,...,K}\G}). This ensures that w, (k=1,2,...K) can have 
roughly the same magnitude as w,U, = wjU; = wjUj. 

For example, assume that the deviational variables dọ and n; in 
the MSCP-MSMS model can vary in different ranges as 
di, ni = [0,120 x 10°], d2, n; = (0, 100 x 109], ds, n3 = [0, 19 x 106], 
da, n4 = [0, 90], ds, ns = [0, 350000], dg, ng = [0, 4000], and d;, n7 = 
[0, 50], resulting in unintentional bias toward objectives. According 
to the above idea, the values of U, and V, can be computed as 
follows: 

U1 4 V4 2 120 x 10°, U; = V4 — 100 x 106, U4 = V4 — 19 x 10°, 
U4 V4 90, Us Vs 350000, Ug = Vg —4000, and U; = 
V;—50. Dy = max (Uk, Vj) — 120 x 105, I; = min (Up, Vi} = 50, 
Ww; = (Lj/Lj-- Dy) = 0.41666 x 1076, wz = (Dy/Lj-- Dy) = 0.999999, 
Ww = (u7U7/U2) = 0.499999 x 107°; W3 = (u7U7/U3) = 0.2631576 
x107, W4 = (u7U7/U4) — 0.555555, Ws = (u7U7/Us) = 0.142857 
x107, and wg = (u;U; /Ug) = 0.0124999. 

This ensures that dx and n, (k = 1,2,...,7) can have roughly the 
same magnitude as w4U, = ... = w;U; ~ 50 for deviational vari- 
ables d, and n, in the objective function of the MSCP-MSMS 
model. The normalized weights applied to the MCGP-MSMS 
model can be expressed as follows: 

(MCGP-MSMS-normalized) 

Min wı(d; +e; )2- wo(dj +e})+w3(d3 +e; )+wa(dy +e, )+ 
ws(ds --ez )-- We(dg --eg )+w7d7 

s.t. Eqs. (3.1)-(3.36). 

The MCGP-MSMS-normalized model is again solved using 
LINGO [28] to obtain the same solution as in the MCGP-MSMS 


Alternatives 


Fig. 5. The solution of MCGP-MSMS-normalized with job creation. 
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Alternatives 


Fig. 6. The solution of MCGP-MSMS-normalized with social acceptance. 


model. In order to demonstrate the usefulness of the proposed 
normalized mechanism, consider the following cases: 

Case I: owing to the increased unemployment rate in Spain, 
DMs hope to increase the employment rate, so the weight added 
to the goal of job creation is increased by 10 times in the MCGP- 
MSMS-normalized model. Solving this problem with LINGO [28] 
yields a new solution, as shown in Table 6 and Fig. 5. As can be 
seen in Table 6, Goals G1, G4, and G5 are achieved, while goals G2, 
G3, G6 and G7 are not in both RC and MCGP-MSMS-normalized 
models. Notably, a value of 72 for the goal of job creation is 
achieved, while the highest achievements of other goals decline 
slightly. It is clear that using the normalized mechanism, DMs can 
define the appropriate weights for the MCGP-MSMS-normalized 
model to achieve their particular goals. 

Case II: due to public resistance, DMs would like to increase the 
social acceptance rate to the maximum possible value, and to do so 
the weight added to the goal of social acceptance is increased by 10 
times in the MCGP-MSMS-normalized model. This problem is 
solved with LINGO [28] to obtain the solution shown in Table 6 
and Fig. 6. Fig. 6 shows the locations for each plant. As seen in 
Table 6, Goals G1, G2, and G5 are achieved, while goals G3, G4, G6 
and G7 are not in both RC and MCGP-MSMS-normalized models. 
Notably, a value of 43 for the goal of social acceptance rate is 
achieved, which is the highest social acceptance rate in this case, 
while the highest achievements of other goals are slightly lower. 
Compared with the RC model, the power generated is 2.62% higher, 
investment cost is 6.97% lower, emissions avoided are 1.22% higher, 
and operation and maintenance costs are 4.53% lower in the MCGP- 
MSMS-normalized model. It is seen that by adding weights to the 
objective function, the proposed method can easily be used as a 
decision aid to determine the best or most appropriate solution to 
multiple objective problems. Each goal of the multiple objective 
problem can be divided to multiple aspiration levels to better suit 
management requirements, such as "the more the better “ or "the 
less the better". These constraints can also be easily added in the 
MCGP-MSMS-normalized model to mirror real-world situations. 
The MCGP-MSMS-normalized model provides a feasible and robust 
way to choose an optimal location for renewable energy plants. 


4. Social acceptance of a wind turbine installation 


Due to the abundant wind resources on the west coast of Taiwan, 
28 possible on-shore locations are evaluated for implementation of 


wind turbine generation systems, as shown in Fig. 7. Though there 
were 314 wind turbine systems installed in 2012, generating 621 MW, 
due to the greater awareness of environmental protection issues and 
participatory democracy, protests against proposed wind turbine 
installation have occurred, such as the public resistance from Yuan Li 
township, Penghu county, XinWu township, Tongxiao township, 
Sanzhi district, Zhunan township, and others. As shown in Table 2, 
there are 886 wind turbines scheduled to be installed by 2030, though 
due to popular protests this will be the most difficult part of the 
subsequent installation of wind turbines in Taiwan. In the past decade, 
there have been many studies on quantitative goals (e.g, maximize 
power generated) in determining the optimal location of renewable 
energy facilities. However, qualitative goals, such as individual pre- 
ferences (e.g, minimize popular protests) have rarely been addressed 
in solutions to this problem. The second objective of this study is to 
provide a feasible method for dealing with the qualitative problem of 
reducing public resistance, in order to facilitate subsequent installation 
of wind power system. For simplicity, the linear utility function is 
given as the following two membership functions (MF), 7;(f;(X)), to 
represent social acceptance rate for the wind turbines installation. A 
graph representation of MF is shown in Fig. 8. (41) is used to 
maximize f;(x)- g; in order to increase the social acceptance rate 
(e.g, the greater the financial compensation, the higher the social 
acceptance rate). 


L d fi(928 

Le if li; «f;(X) «gi, (the right linear MF) (4.1) 
0 if | fix <li 

(4.2) is used to minimize f;(x) £ g, in order to increase the social 


acceptance rate (e.g., the less the wind turbine noise, the higher the 
social acceptance rate). 


nifi) = 


1L, f fis 
nif (xX) = ae if gf < ui, (the left linear MF) (4.2) 
0o if fizu 


where u; and l; are the respective upper and lower tolerance limits for 
the ith goal. The achievement model with the MFs, can be formulated 
as follows: 

(MCGP-U) 


Min Y laidi +d; )2-B(ej +e-)+4f; +wh; ] 
ici 


s.t. (fix) Sa the right linear MF 


fi(x)—d;* +d; =g;, the upper limit value g; ẹ for fix) > gi 
nif (x)) +f; — 1, the highest value of 1 
nif) = 4-2, for left linear MF 
fico ej er =g;, the lower limit value g; for fi(x) < gy 


nf (%)+h; — 1, the highest value of 1 


dt, dy, fj, ufi) = 0 


xeF (F is a feasible set) 
where @;,f;,4;, and y; are the weights attached to the deviational 
variables, d; d; ‚et e; f; , and hj . 

In general, Taiwan citizens are opposed to high-polluting coal- 
fired and fossil fuel power generation, and the safety of nuclear 
power is also an issue of public concern, so there is support for 
renewable energy. But residents still do not want to have wind 
turbines sited too closed to their dwellings, schools or workplaces 
since wind turbines are a source of noise in rural environments. 
The most common complaints and source of social resistance to 
wind farms are noise and blade flicker (glare reflecting off the 
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spinning turbine blades). Especially, the lower frequencies of noise 
from the wind turbine are judged to be very annoying during the 
night. These two issues are directly related to the distance from 
the wind turbines to residents homes. Thus people would like to 
maximize the distance between wind turbine installations and 
their dwelling and workplaces, to achieve higher satisfaction. 

But due to the costs associated with construction and land 
acquisition, and considering Taiwan's dense population, it is often 
difficult to fully meet the requirements of residents. So, finding a 
balance between the requirements of residents and construction 
considerations is the most important issue. For sake of simplicity, 
three basic attributes are taken into consideration: power gener- 
ated (MW), investment cost (I), and social acceptance (S) with MF. 
Table 8 shows the data for MW and I attributes. Table 9 shows five 
candidate locations for the installation of wind turbines. Table 10 
shows the social acceptance with MF. Based on (4.1) and Table 10, 


Table 8 
Alternatives for electric generation. 


Alternative Variables MW I (US$) 
Wind turbine 1 (1 <P < 2 MW) x 1.98 3,280,000 
Wind turbine 2 (2 <P <4 MW) X 3.92 6,560,000 
Wind turbine 3 (4 < P < 10 MW) x3 9.12 9,840,000 
Wind turbine 4 (10 < P < 15 MW) X4 14.55 10,500,000 
Table 9 
Five places attributes. 
Places Distance to residents house (m) Land cost (US$) 
5 300 70,000 
6 500 90,000 
7 600 98,000 
8 1200 140,000 
9 450 80,000 
Table 10 


Wind turbine (noise vs. social acceptance). 
Alternative Distance noise 


Wind turbine 1 (1 < P <2 MW) [200 m,1000 m]=[40 db,30 db] 


Wind turbine 2 (2<P<4 MW) 


[250 m,1200 m] =[42 db,35 db] 


Wind turbine 3 (4<P<10MW) 


[300 m,1500 m] =[45 db,38 db] 


Wind turbine 4 (10 < P < 15 MW) [400 m,1500 m] =[46 db,39 db] 


[0,1] (4.1) right linear MF 


[0,1] (4.1) right linear MF 


[0,1] (4.1) right linear MF 


[0,1] (4.1) right linear MF 


all possible levels of social acceptance (satisfaction) for wind 
turbine locations are calculated, as shown in Table 11. Four 
possible wind turbines installations (see Table 8) for power 
generation in five possible locations (see Table 9) must be 
evaluated in order to satisfy the following goals. 


- (G1) Power generated must be higher than 15 MW and the 
higher the better. 

- (G2) Investment cost must be limited to 50 x 10°$/year and the 
lower the better. 

- (G3) Social acceptance must be limited to the level of 50 and 
the higher the better. 


Table 11 
Satisfaction of wind turbine installation. 


Alternative Place Value of MF (satisfaction) 
Wind turbine 1 5 0.125 
6 0.375 
7 0.5 
8 1 
9 0.3125 
Wind turbine 2 5 0.052 
6 0.263 
7 0.368 
8 1 
9 0.21 
Wind turbine 3 5 0 
6 0.166 
7 0.25 
8 0.75 
9 0.125 
Wind turbine 4 5 0 
6 0.09 
7 0.18 
8 0.727 
9 0.045 
Social acceptance Corresponding MF Graphs 


Satisfaction 


E 


Distance 
200m 1000m 


Satisfaction 


L^. 
0 Distance 


250m 1200m 
Satisfaction 


0 Distance 


300m 1500m 
Satisfaction 


0 Distance 


400m 1500m 


C.-T. Chang / Renewable and Sustainable Energy Reviews 41 (2015) 379-389 387 


Alternatives Places 


8 
QO OW 
c2 O00 

Q 


Fig. 9. The solution of MCGP-U without weight. 


Alternatives 


A 
0 


BS 


Fig. 10. The solution of MCGP-U with weight. 


The achievement model without adding weights to each goal is 
formulated using the proposed model MCGP-U (see the Appendix 
for details). The model is then solved using LINGO [28] to obtain 
the solution as (MW, I, S)=(29.57, 34.26 x 10°, 0.4287). Fig. 9 
shows the locations for each wind turbine installation. As the 
result the following goals are achieved: social acceptance of 
42.87%, power generation of 197%, investment cost of 131.5% 
achieved. Clearly, the problem can easily be solved using the 
proposed MCGP-U. In order to increase the social acceptance rate, 
a large weight is added to G3. Then, the problem is solved again to 
obtain solution as (MW, I, S)=(16.53, 16.16 x 10°, 0.6135). Fig. 10 
shows the places where each wind turbine should be located for 
higher social acceptability. As seen in this solution, the social 
acceptance rate increases from 42.87% to 61.35% when the large 
weight is increased to G3. This shows that the proposed model 
provides feasible features for a DM to deal with multiple decision 
making problems. 


5. Conclusions 


Mitigation of climate change and environmental pollution has 
increased the need to replace fossil fuel energy with renewable 


energy world-wide. The capacity expansion planning problem of 
renewable energy industry is an important issue for both devel- 
oping and developed counties in the implementation of appro- 
priate renewable energy policies. To deal with this problem, 
Ramón and Cristóbal [9] proposed a method in which an optimal 
outcome of plant types, numbers of plants and locations can be 
easily determined. In this study, a MCGP model is proposed to 
improve the usefulness of that model of through an example of 
locating five renewable energy plants for power generation in five 
locations in the autonomous region of Cantabria, in the north of 
Spain. Moreover, the normalization method is introduced to the 
MCGP for helping DMs to set appropriate weights according to 
their preferences or policies. As a result, MCGP can be used as an 
efficient, effective and feasible tool for DMs when dealing with the 
capacity expansion planning problem of renewable energy indus- 
try. In addition, social acceptance of a wind turbine installation in 
Taiwan is also considered. 


Appendix 


(MCGP-MSMS) 

Mind; +e +d; +e}+d; +e; + dy +e; +d; +eł + dg + 
eg +d} 

s.t.—The following for power generated goal (the case of the 
more the better): 


60, 959, OOOF; 8 +... +37, 770,000 F;14—dj +d; =y}, (3.1) 
yi-et +e; - 120 x 106, (3.2) 


110 x 106 < y, < 120 x 10°, the sufficiently large value for power 


generated is 120 x 10° W (3.3) 
The following for investment (the case of the less the better): 
23,425, 000F, s +... +47, 958, 000F7 14— dj +d; — y», (3.4) 
y;—ej +e; —80 x 10°, (3.5) 


80 x 10° < y, < 100 x 10°, the sufficiently small value for investment is 80 x 10° 


(3.6) 


The following for emission avoided goal (the case of the more 
the better): 


9, 649, 680F1 8 + ... +4, 839, 548F7 14 — d7 Hd; — ys, (3.7) 
ys—ej +e; —19 x 10°, (3.8) 


18 x 10° < y, < 19 x 10°, the sufficiently large value for emission 
avoided is 19 x 10t (3.9) 


The following for jobs created goal (the case of the more the 
better): 


15F18 +... +20F714- di +d} = y4, (3.10) 


y4—e4 +e, =90, (3.11) 


70 x y4 x 90, the sufficiently large value for jobs created is 90 
(3.12) 


The following for operation and maintenance costs goal (the 
case of the less the better): 


37750F; s +... +75,000F714-d3 +d; —ys, (3.13) 


ys—e¢ +es — 300000, (3.14) 
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300000 < y; < 350000, the sufficiently small value for OM costs is 300000 


(3.15) 
The following for total distance between nodes goal (the case of 
the more the better): 


6.5FasiFosi t+... + 75F 135iF 14si — dg; +dg = yg. (3.16) 


Yg eg +e; = 4000, (3.17) 


2910 x yg x 4000, the sufficiently large value for total distance is 4000 
(3.18) 


The following for social acceptance goal (the case of the more 
the better): 


7Fig+...+9F714—d7 +d; =50, (3.19) 
Fsosi — Fso1 — Fso2 — Fso3 — Fso4 — Fsos — Fso6 — Fso7 = 0, (3.20) 
Fgsi + Fosi + Frosi + Fi isit Fiasit F13si +F14si — Fsiso = 9, (3.21) 
Fig +Fog + F3g + Fag — Fesi = 0, (3.22) 
F29 + F49 + Fgo — Fosi = 0, (3.23) 
F110 +F310 + F510 + Feo — Fiosi = 0, (3.24) 
F111 +F311 +F511 Fei +F711 —Fiisi = 0, (3.25) 
F212  Fsi2  F712 — Frasi = 0, (3.26) 
F213 + Fa13 +F613 — Fiji = 0, (3.27) 
Fyya+F314+F14—Fyasi = 0, (3.28) 
Fso1 — Fig — F110 — F114 = 0, (3.29) 
F502 — Fag — F29 — F212 — F213 = 0, (3.30) 
F503 — F38 — F310 — F311 — F314 = 0, (3.31) 
Fso4 —F4g — Fag — F413 = 0, (3.32) 
Fsos — F510 — F511 — F512 = 0, (3.33) 
Fso6 — Feo — F610 — F611 — F613 = 0, (3.34) 
Fso; — F714 — F712 — F714 = 0, (3.35) 
F gsi + Fosi + Prosi + Frisi + Fiasi + Fi3si + Frasi = 5, (3.36) 


where Eqs. (3.1)-(3.19) correspond to the goals power, investment, 
emissions, jobs, operation and maintenance, and social accep- 
tance; and Eqs. (3.20)-(3.36) are the conservation of flow con- 
straints for each node. 

(MCGP-U) 


Mind; +d} +e; +e} +d; +d} +e; He) 


s.t. —The following for power generated goal (the case of the more 
the better): 


1.98F15+1.98Fi6+...+14.55Fs9—dj +d; =y; (4.1) 
10 < y, < 30, (4.2) 
yi-et +e; =30, (4.3) 


The sufficiently large value for power generated is 30 MW . (4.4) 


The following for investment (the case of the less the better): 


(3.280, 000 +700, 000)F 15 +... + (10, 500, 000 + 800, 000)F49 -d7 --dj —y», 


(4.5) 
5,000, 000 < y; < 10,000, 000, (4.6) 
y;-—ej re; = 5,000, 000, (4.7) 


The following for social acceptance (the case of the more the 
better): 


0.125F15 +0.375Fi¢+...+0.045F 49 


di +d; = 48 

Fi5+Fig+...+Fag dii isnt Oel 

0.5 £ y3 <1, (4.9) 
y3—e3 +63 =1, (4.10) 


The following for one wind turbine can only assign one place, 
and vice versa: 


9 

£ Fy<1, i=1,2,3,4 (4.11) 
j=5 

4 

È Fj<1, f=1,2,3,4,5 (4.12) 
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